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In Euclidean four-dimensional S\J{N) pure gauge theory, eigenvalue distributions of Wilson loop 
parallel transport matrices around closed spacetime curves show non-analytic behavior (a 'large- 
phase transition') at a critical size of the curve. We focus mainly on an observable composed 
of traces of the Wilson loop operator in all totally antisymmetric representations, which is regu- 
larized with the help of smearing. By studying sequences of square Wilson loops on a hypercubic 
lattice with standard Wilson action, it is shown that this observable has a nontrivial continuum 
limit as a function of the physical size of the loop. We furthermore present (preliminary) numer- 
ical results confirming that, for large A^, the dependence in the critical regime is governed by 
the universal exponents 1/2 and 3/4 as expected (Burgers universality). 
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1. Introduction 

Already in 1981, Durhuus and Olesen [1] discovered that Wilson loops in pure SU(A'^) gauge 
theory in two Euclidean dimensions exhibit a transition from an 'ordered phase' to a 'disordered 
phase' at large N. In two dimensions, the eigenvalue spectrum of the untraced Wilson loop unitary 
matrix depends only on the area that is enclosed by the (nonselfintersecting) space-time curve 
defining the loop and a sharp infinite-A^ phase transition occurs at a critical value of this area. At 
this point, the gap in the spectrum, that is present for small loops, closes in a non-analytic way. 

In Ref. [2], Narayanan and Neuberger conjectured that a similar large -A/^ transition also occurs 
in three and four dimensions and that all transitions belong to a single universality class. If this is 
indeed the case, it might be possible to connect the perturbative regime with non-perturbative mod- 
els in four dimensions using the universal properties of the transition at a critical scale. However, 
so far the universality conjecture has been confirmed only for the three-dimensional case [3]. 

In the following, we present a numerical study of the transition and its universal properties in 
four dimensions using lattice methods. 

2. Main observable 

The Wilson loop matrix associated with a closed spacetime curve is defined as usual: 



where all totally antisymmetric representations of W enter (/ denotes the fundamental representa- 
tion and Xk"^"^ denotes the character in the ^-fold totally antisymmetric representation of SU(A'^)). 
The variable y is defined such that a zero of the determinant at j = corresponds to an eigenvalue- 
angle Q = ±71 on the complex unit circle. The gap in the eigenvalue spectrum closes at = ±;r 
(cf. Fig. 1 below). Therefore, useful information about the large-A'^ non-analyticity can be obtained 
by expanding ^^(j, in powers of y around y = (with Qcp = 0, only even powers of y enter): 




(2.1) 



with path-ordering operator I^. 

To study the large-A^ transition in the spectrum of W , we focus on 



^^(3;,^) = (det(.« +e-mf{^))) = Ie("-^)>'(zr™(W(^))> 



k=Q 



(2.2) 



(j, ^) = ao (^) + a 1 C^)3;2 + ^2 (^)/ + O ( j^) . 



(2.3) 



In particular, it turns out (see below) that the scaling-invariant ratio 



aoC^)a2('r) 



(2.4) 



provides a valuable signal for the large-A^ phase transition in the spectrum of W . 
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3. Universality class from heat-kernel model 

The heat-kernel probability density (w.r.t. the Haar measure) for an SU(A'^) matrix W is given 



by 



implying 



all irred. r 



(3.1) 



(3.2) 



due to character orthogonality. The parameter t can be interpreted as a diffusion time, dy and C2(r) 
denote the dimension and the quadratic Casimir invariant of the irreducible representation r. The 
above probability density can be realized in a simple multiplicative random matrix model and holds 
exactly in two-dimensional pure Yang-Mills theory when t is identified with the dimensionless area 
variable g^Ns^ (with 't Hooft coupling g^N and denoting the area enclosed by the loop). The 
heat-kernel single eigenvalue-angle distribution is given by [4] 



l^{-\Yl^d{p,q)e-^'^^P''^co^{{p + q+l)B) 

{N+q)\ 



(3.3) 



p=0 ^=0 

withC(p,^) = \{p + q+\) (N-l^+q-p) andd{p,q) - p,g,^N-p-iy. p+q+i 



An immediate consequence of the group theoretical structure (3.1) is that 



■log<^(3;,T) 



log( det (e2 +e '^-Wf 



HK 



(3.4) 



Ndy ^ ' Ndy 

with T = f (1 + 1 /N) satisfies Burgers' equation [5] 

with initial condition (l)^^{y,0) = — ^ tanh ^. At A'^ = oo, Burgers' equation produces a 'shock- wave' 
singularity at 3^ = when T reaches the critical value Tc = 4, cf. Fig. 1. 
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Figure 1: (left) and heat-kernel single eigenvalue density (right) in the limit N ^ °° for various t 



From Burgers' equation (at finite A'^), we get 



(3.5) 
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At infinite N, the inverse slope of at j = increases from —4 at T = to at T = 4 since 



1/2, 0<T<4, 
1/6, T>4, 



resulting in a discontinuous jump in for T > 4 at y = 0. This discontinuity results in a 

nonzero single eigenvalue density at = ±71 (on the unit circle in the complex plane). The infinite- 
N phase transition (non-analytic behavior with chai^acteristic exponents) in the eigenvalue density 
of W occurs at the point where the gap in the spectrum disappears (this happens at T = = 4), 
cf . Fig. 1 . 
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Figure 2: (0^^{t) for various A^. 



The singularity is absent for any finite A'^, and very large values of A'^ are needed to see its 
formation, cf. Fig. 2. From Burgers' equation we obtain: 



3 a\ 

lim N^i — 

N^oo ClQ 



3 a? 

lim A^"2 — 



ai 



T=4 



T=4 



lim (ON\r-A 



1 3 1 

24 V 2 
1 



K: 



At? 



1.046, 



lim N 2 — — 



1 



T=4 



6 V 2 



K{K^-\) 



(3.6) 

(3.7) 
(3.8) 
(3.9) 



Furthermore, the (purely imaginary) roots of G"]^ {y, t) in the critical regime (around j = 0, T = 4) 
scale like ? . 



4. Numerical results in 4D 



In four dimensions, the Wilson loop operator develops a perimeter divergence (and additional 
corner divergences if the spacetime curve has kinks). A convenient way of regularization is by 
smearing (cf. Ref. [6] for details). We deal with the intrinsic UV divergences of the action by using 
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Figure 3: Lattice results for (Oig{b,L) (with cubic spline interpolation between data points). 
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Figure 4: mib^L) and TN{b,L) foxN^\9 (solid) andA^= 11 (dashed) and L = 2,4,6, 8. 



a lattice discretization of the theory. To extrapolate to the continuum, we compute sequences of 
square Wilson loops of sides 1 < L < 9 for inverse 't Hooft couplings 0.348 <b = < 0.374 at 
A'^ = 19, using a single plaquette Wilson action on a hypercubic lattice (volumes 12^ and 14^ are 
used to exclude data contaminated by finite-volume effects) and a combination of heat-bath and 
overrelaxation updates. The amount of smearing S is chosen proportional to the size of the loop, 
5 = L^/55. We then measure C0N{b,L) (determined from the average characteristic polynomial) as 
a signal for the large -A/^ phase transition. Laige,-N scaling is tested with additional runs a.tN = 11. 

Figure 3 shows our numerical results for (0\'){b,L). Similar to the heat-kernel model, we 
observe that very large values of A'^ would be needed to allow for a direct observation of the singular 
large-A'^ behavior in (Om- Therefore, we define an A'^-dependent map (ON{b,L) — ;> ^^{b.L) through 

w^(ft,L) = w^'^(T^(^,L)). (4.1) 

The required inversion is unique and the map can be used even though the r-dependence of {Xr{W)) 
in 4D differs from exact Casimir scaling. While the formation of the jump in (Oj^ is slow, 
converges rapidly to x^, cf. Fig. 4. This implies that, to some subleading order in ^, (ON{b,L) « 
w^^{r^{b,L)). 

This relation can be taken over to the continuum limit, where the two variables b and L get 
replaced by a single length variable / = L/Lc{b), the side of the loop in physical units. To set the 
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Figure 5: Example for continuum extrapolation at = 19 and log(^) = —0.2. 
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Figure 6: Continuum functions %(/)^^ for = 19 (black) and = 1 1 (brown). The gray bands show the 
corresponding error estimates (jackknife), the horizontal dashed line corresponds to the critical value Tc = 4. 

scale, we use the critical deconfinement temperature 1 /Lc{b), determined by 

The continuum limit of % is obtained by extrapolating XN{b,L) at fixed / to L — )• oo, cf. Fig. 5. 

We observe that % has a nontrivial continuum limit which is a smooth function of the physical 
loop size / for all A'^, cf. Fig. 6. This establishes the transition and its universality since we can re- 
place (^A?(y, ^) by i^^^{y^ "^"(O) the vicinity of the critical point (y = 0, T = 4) without changing 
the singular large-A'^ properties. The dependence on / is consistent with asymptotic freedom as T 
modulo a shape-dependent factor can be interpreted as an effective coupling constant. 

Additional checks of universality can be obtained from the critical exponent 3/2 of the ratio 
ao/a\ and the A'^^^/^-scaling of the level density in the critical region (cf. Fig. 7). Both checks work 
very well: numerically, we obtain exponents of 1.52 and 0.73, respectively. 



5. Conclusions 

We have obtained numerical evidence, by extrapolating results of lattice simulations to the 



6 



Large-N phase transition of smeared Wilson loops in 4D pure YM 



Robert Lohmayer 




1/N I 

0.00 0.02 0.04 0.06 0.08 0.10 3 5 6 7 8 



Figure 7: Left: appropriately scaled angular difference A0 between the two peaks in the single eigenvalue 
density closest to = ;r at T = 4 as a function of (blue: data points, red: heat-kernel value). 
Right: N^l^^e as a function of T (blue-green points for = 19, red-yellow points for = 11; black lines 
show the corresponding heat-kernel functions). At the intersection point (at T = 4), A0 ^ N^^^^. 

continuum, that smeared Wilson loops in 4D continuum pure SU(A'^) gauge theory undergo a large- 
N phase transition at a critical loop size. Furthermore, we have confirmed the expected universal 
properties of the transition. For a complete presentation of our results (including data for larger N) 
we refer to Refs. [7, 8]. 
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